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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

APPLICATION OF STATISTICAL THEORY TO BEAM-RIDER 
GUIDANCE IN- TEE PRESENCE CP NOISE. 

I - WIENER FILTER THEORY 
By Elwood C. S-tewaorb 

SUMMARY 


A study has been made of the application of Wiener filter theory to 
the design of a beam-rider guidance systea operating in the presence of 
glint noise. Target and missile motions are restricted to the same plane. 
The Wiener theory is then used to establish the theoretical lo-wer limit 
of root-mean-square error and the coire spending desired transfer- function 
charac-fceristicB . It is shown that although the practical achievement of 
these results is restricted by saturation effects, the theory is useful, 
with suitable modifications, as a guide in system design. Such modifica- 
tions have been applied to the design of systems for which the optimum 
filtering is placed either in Hhe miss lie- control system, or in the tracking 
radar. The error performance of -these sys-tems for different noise magni- 
tudes is presented. Other considerations such as servo energy requirements 
are briefly discxissed. 


INTRODUCTION 


In the design of a missile-guidance system certain standard criteria, 
such as the fastest response, are not al-ways the most useful. This is 
particularly true in the case of guidance systems -which are forced to 
operate in the presence of certain random imwanted disturbances known as 
noise. The reason for -this is that the effects of the noise can seriously 
reduce the probability that the missile will hit the target. Furthermore, 
■the sources of noise, being dependent on the target characteristics, can- 
not be eliminated. Consequently the guidance systen shovild be designed 
to minimize the miss distance even -when the noise is present. This prob- 
lem will be considered here. 

There are -two possible design approaches. In the first, the form 
of the system’ (i.e. , the transfer function)- is assumed and an attempt is 
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made to adjust the existing parameters so as to reduce the effects of 
noise (ref. l) . This procediare is not only difficult to apply, but the 
ultimate performance is limited by the assumed form. In the second 
approach, the form as well as the parameters of the transfer function 
are determined so that the noise has the' least possible effect on the 
performance of the missile; an attempt is then made to devise a system 
which has a transfer f;mction approximating this optimum transfer func- 
tion. The latter approach will be adopted here. 

The problem of determining optimum transfer functions has been 
encountered previously in the communication field, and as a result of 
this encounter a statistical theory known as Wiener filter theory has 
been developed. By the use of this theory it is possible to determine a 
unique optimum transfer function which -will result in a theoretical, lower 
limit of mean-sq\aare error between the desired and the actual missile 
position. Very little work has been done: in the application of this 
theory to beam- rider guidance. Previous works, references 2 and 3> have 
been confined to simple homing systems. The pui^ose of the present paper, 
however, will be to Investigate the applicability of this theory to a 
beam- rider guidance system. ‘ 

In the application of this theory to missile guidance it is necessary 
to make certain assumptions. Foremost of . these is the assumption that the 
target and missile move in the same plane, taken in this report to be 
horizontal. Other assumptions, such as those relating to the class of 
target maneuvers and noise, are discussed in the text. Within these 
restrictions, however, the theory may be used to obtain a measure both 
of the error performance that might be expected and the difficulties to 
be overcome in order to realize this performance. 


SYMBOLS 

R noise magnitude or zero frequency spectral density, ft^/r adlan/sec 

Tjj time constant of the noise spectrum shaping filter, sec 

Yq optimum closed-loop transfer fiinction_. 

a acceleration of target maneuver, ft/sec^ 

k twice the average switching rate of target acceleration, l/sec 

y apparent target displacement from true target center due to noise, 

^ ft 

y^j missile displacement from a space reference, ft 
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•barget dlsplacsnent frcm. a space reference, ft 

spectral density of noise displaeeinent ft^/radian/sec 

spectral density of target displacement y^, ft^/radlan/sec 

error between target and missile position, y^ - yj^, ft 

component of error € due to noise, ft 

ccmponent of error e due to target motion, ft 

control-surface deflection, radians 

angle of yaw, radians 

optimum open- loop transfer function 

open- loop transfer function of system approximation to Pq 
angular frequency, radians/sec 


GENERAL CONSIDERATIONB 


Of the many sources of noise tdiich may exist in a guidance system 
utilizing radar detection, glint noise is one of the most serious. Glint 
noise is a term that is used to describe a shift in the apparent center 
of a target as determined by a tracking radar. It is due basically to 
the variable reflection charaeteristics of aircraft targets and arises 
from the relative movement of the various reflecting surfaces. Since the 
radar utilizes the reflected signal to determine target position, vari- 
ations in the reflected signal are interpreted by the radar as a shift 
in the target center. This type of noise is . parti ciolarly Important since 
it is due fundamentally to the target characteristics and therefore cannot 
be eliminated by any known radar improvements. The sittxatlon is illus- 
trated in figxrre l(a) where the true target position is Indicated as y^, 
and the glint noise is represented by the displacenent yjj. 

The present ^report is restricted, for the sake of simplicity, to a 
two-dimensional study in which the target and missile move in a horizontal 
plane. ^ The guidance system is considered to be of the beam-rider t 3 ?pe, 
illustrated in figures l(a) and l(b) . It should be noted here that dis- 
placements are referred to a fixed apace reference. The function of the 
guidance system is to make the missile position yj^ coincide as closely 

^The cfmplet^. three- d-tme na-i nnal p-rnblem wnnl d -requii^ a mn-re crgnpl^x 
analysis than used herein. Possibly either the present theory or Wiener’s 
theory for multiple time series Cref . 4) could be applied to this case. 
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as possible with the actizal target position This requirement would 

be relatively easy to satisfy if the tracking radar could locate the 
target precisely. However, because of the reidar noise the only informa- 
tion available to the guidance syst^n is the apparent target position as 
illustrated in figui^ iCb) . For this reason the task is much more diffi- 
cult and the missile position may deviate considerably from the actual 
target position. The dlffei^nce, ym - yw, is denoted on this figure by 
the error e, which should be minimized in some manner. 

A suitable criterion for Judging system performance depends primarily 
on the manner in which the system operates. In the case of the beam-rider 
system the missile -to- target range is not normally transmitted to the 
missile so that the missile cannot know when the target will be reached. 
Hence the error should be minimized for all values of range or, equiva- 
lently, of time. A mathematically convenient criterion which does not 
involve weighting with respect to missile travel time is the mean-square 
time average of the error. This criterion will be used herein. 

The design of a system normally depends on the inputs to be encoim- 
tered, in this case the target motion and noise. Because of their randcan 
nature it is not convenient to define these quantities explicitly aa 
functions of time, and statistical descriptions ere more suitable. Since 
it is generally believed that the target motion and noise are uncorrelated 
they will be described independently. What follows is a brief discussion 
of these inputs . 

Intensive effort has been devoted in recent years to the measurement 
of radar glint noise. References 3 5 tiirough 8 are typical of such 

work. The quantity of most general Interest in these measxirements is the 
displacement of the apparent center of the target from the true center, 
or shown in figure 1(a). This quantity can be defined statistically 

by means of (l) the amplitude distribution and ( 2 ) the power spectral den- 
sity. Altho\ 2 gh the determination of these quantities is somewhat uncertain, 
it is generally found that the amplitude distribution is approximately 
Gaussian and that the spectral density can be adequately represented by 


= 


N 


( 1 ) 


Spectra obtained from any one individual experiment may deviate somewhat 
frcaa this form but it is generally considered that the above character- 
istic represents a reasonable average of many different experiments. 
Examination of noise spectra Indicates that the break point (l/2«Tjj) is 
on the order of several cycles per second and for this report will be 
taken to be 6 cps corresponding to Tjj - O.O 265 second. The magnitude 
of the spectrum, N, depends on factors such as target size and target 
aspect so that the gtiidance system is usually forced to operate over a 
wide range of magnitudes . This range may extend frcmi 7 ft^/radian/sec 
for small targets up to around 30 ft^/radian/sec for large bombers. 



MCA EM A55E11 


5 


Becaxise of the variation in magnittidg due to uncertainty in target size 
and aspect as well as uncertainty in the noise measurements, this factor 
becomes of real and practical importance and will he considered herein. 

In considering target maneuvers it is difficult to say exactly how 
a target will maneuver when under attack. However, a reasonable situa- 
tion might be one in which the target is merely aware of the attack and 
therefore maneuvers in some severe manner to avoid being hit. Here it 
will be assimed that the target maneuvers laterally with maximum acceler- 
ation alternately in opposite directions. The duration of each acceler- 
ation will be a randcsa function determined by some distribution. A 
reasonable distribution which leads to an easily handled spectral density 
is the Poisson distribution [ (l/T)exp(-T/T) ] where T represents the time 
and T the average time between switching of the acceleration. As shown 
in reference 9 "the spectral density of the target acceleration is then 
described by 


s= ka^ 

+ k^) 

Here the quantity a represents the magnitude of the target acceleration 
normal to the beam, and k is twice the average switching rate, k = 2/T. 
The spectral density of the target displacement is then given by 

It would appear that there is a problem here concerning the existence of 
this spectral density because of the in the denominator. However, 

it can be shown that it is possible to use hhis representation for pur- 
poses of computations (see Appendix A) . For the tail-chase maneuver to 
be used in a later example the target acceleration is specified to be 
±1 g at an average period of 5 seconds, which gives a = 32.2 ft/sec^ 
and k = 0.4 switch/sec. 

It should be pointed out that a system design based on the target 
motions described above would operate well against this class of maneuver 
as a whole. This appears to be a desirable procediare. Nevertheless, 
without altering the parameters this system would not be expected to 
operate as well as it could against one particular target maneuver such 
as a single target turn. Even then, however, it can be shown by simu- 
lation studies that systems optimized for the statistical maneuver used 
herein are essentially optimum for the single turn maneuver as well. 
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ANALYSIS, RESULTS, AND DISCUSSION 
The Wiener Filter Theory 


The prohlem of minimizing the effects of noise can be considered to 
be one of compromise. At one extreme for which the system response is 
fast the error becomes large because of the ability of the missile to 
follow the noise too well. At the other extreme for which the system 
response is too slow the error becomes excessive because of the difficulty 
the missile has in following the target maneuvers . The optimum system is 
one which will compromise these two situations in the best possible man- 
ner, More precisely stated the problem beccxnes: Given the statistical 

characteristics of the two input quantities, target motion and noise, 
what is the optimum .transfer function which will minimize the mean-square 
error The answer to this problem can be determined by a statistical 

theory known as the Wiener filter theory. The final result of this theory 
is an integral which when evaliiated represents the optimum linear transfer 
function. This transfer function is given by 


YoCit*j) 


1 




®"(a) 


do, dt 


( 3 ) 


In this equation 4"^ and 4>" have the following meaning: If $ is the 

spectral density defined by the equation 




then and must satisfy 


(j)” = <l> 

where has poles and zeroes only in the upper half of the complex u 

plane and only in the lower. More details are given in Appendix A. 

The derivation of equation ( 3 ) is beyond the scope of this report; for 
this derivation see references 4 and 10. 

The transfer function Yq of equation ( 3 ) is a mathematical repre- 
sentation of the box in figure 1(b). According to the method of deriva- 
tion, the transfer function must be physical 1 y realizable which means 
that the system is not required to respond to an input before that input 
occurs . 

It might be pointed out that a general solution exists which involves 
the problem of prediction as well as filtering.,. However, since it is. 
apparent from figure l(a) that we are interested only in the missile 
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position coinciding with the present target position, the prediction aspect 
has been eliminated from the general solution in writing equation (3) . 

Bestrictions involved in the theojry .- There are certain restrictions 
implicit in the derivation of equation (3) bo that the validity of its 
application to the heam-rider guidance system depends on how well these 
restrictions are met. First, the input quantities, target motion and noise, 
must be stationary randcmi series (see ref. 1 for a detailed definition) and 
he defined by corresponding power spectra which are continuous. A discon- 
tinuity in the spectnmi might be due to a predictable ccanponent such as a 
sine wave; such components must be eliminated from the input before the 
theory can be applied. It is generally believed that displacements at the 
target are approximately stationary random series (ref. 3 )* Since the 
beam-rider system operates from these displacements, the Inputs to the 
beam- rider system are also stationary random, series. (In contrast, the 
inputs to a proportional-navigation guidance system are nonstationary ran- 
dom series since angular inputs are measured by the missile itself and the 
angles tend to become larger as the range decreases . ) Second, the solution 
is based on linear theory and furthermore is applicable only to a system 
with constant coefficients in its differential equation. On the other 
hand, the kinematic loop of guidance systems generally involves a time- 
variable range factor. In particular for the beam-rider system, the time- 
variable factor is the ratio of the launcher-to-missile and the launcher- 
to- target ranges. However, since the miss is determined primarily by what 
happens near the end of fli^t, during \diich the variation in this ratio 
is small, it is reasonable to assume that the requirement of constant 
coefficients is approximately met. 


Evaluation of the optimum, transfer fiinction .- Within the above 
restrictions, it is possible to use equation (3) to evalioate first, the 
general form of the optimum transfer function and second, the numerical 
constants. As an illustration of this method the theory has been applied 
with certain simplifying assumptions as discussed in Appendix A to the 
target motion and noise characteristics described earlier. For this case, 
then, equation (3) can be evaluated to give the general form of the optimum 
closed- loop and corresponding open- loop transfer fimctions as follows: 


+ 1 ) 

(TpS + l)(Ty2g2 ^ 2^^TyS + 1) 


(^) 


l^o(s) 


- Hi 


(T^gsg + ^^TgS + 1) 

s2(t^8 + 1) 


( 5 .) 


These equations can be evaluated for any specific case and, as an exan^le, 
for the specific target motion and noise char acteid sties given below 
equation ( 2 ) and with the magnitude of the noise, N, chosen to be 
15 ft^/radi an/sec, the optimum transfer functions become 
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(0.8^^ + lAl s + 1) 

(0.687 s + 1) (0.1^90 s® + 0.727 s + 1) 


( 6 ) 




7 kp ( 0 . 8 ^^ 3^ + 1 .^ 3 - 6 + 1) 
s 2(2,50 b + 1) 


( 7 ) 


The chosen value of N represents a mid-value between the expected 
extremes. The significance of this choice will he discussed later. 


The corresponding frequency characteristics are possibly more llliis- 
trative to the control designer. The solid curves in figure 2 (a) show 
the optimum closed-loop characteristic for the above conditions and 
describe the characteristics of the box in figure l(b) . It might be noted 
that since the break point of the noise spectirum occtirs at 6 cps, it is 
essentially flat over the frequency range of Importeince of the optimum 
transfer function. Since actual systems are usually designed on an open- 
loop basis, the corresponding open-loop characteristic is shown in fig- 
ure 2 (b) , The considerations Involved in achieving the characteristics 
of figure 2 are discussed in later sections. 


RMS error performance .- Although the Wiener filter theory can be used 
to define the optimum transfer function, it does not give the minimum 
error directly. This must be evaluated from the optimum transfer function. 
With the earlier aas-umptlons as to the form of the target motion and noise, 
reference 1 shows that the total mean- square error is composed of target 
and noise components which can be evaluated according to equation (8) . 


e2 = 





Yo(lw) |^<&rp(cj)du) 



Yo(lw) 


( 8 ) 


The quantity Yo(iw) is obtained from equation ( 4 ) by the substitution 
s = iw. It is possible to evaluate these integrals mathematically but 
experience has shown that it is easier and more Instructive to use 
graphical techniques. Tor the exan^le iinder discussion, evaluation of 
equation ( 8 ) gives: 

^/e^ = 13.7 feet 
•J = 6.7 feet 
= 15.^ feet 
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The later figure represents, then, the theoretical lower limit of rms 
error that can he acMeved hy a linear constant-coefficient transfer 
function cori*esponding to the assumed target and noise characteristics. 

For later tise, the integrands of equation (8) for the previous 
example are plotted in figure 3 - Physical 1 y, these curves can he inter- 
preted as the frequency distrihution of the energy in the components of 
error. It is seen that the error energy is concentrated in the frequency 
range helow 2 cps. 

As mentioned earlier, it is important to consider a large range of 
noise magnitudes for a nijmher of reasons. For exaii 5 )le, measurements of 
glint noise are subject to considerahle discrepancy and may cover a wide 
range of magnitudes. Likewise, the noise magnitude may vary due to 
different sized targets or different attack aspects. In order to assess 
these effects the Tn1r>-TTmnn error has heen determined as a function of the 
zero frequency noise magnitude N of equation (l) hy a procedure identical 
to that above. Thi s requires the determination of a different transfer 
function Yq for each value of IT, as in Appendix A for N = 15 ft^/ 
radian/sec, and the evaluation of the resultant error hy equation (8). 

He result is shown in figure 4 hy curve A. Each point on this curve 
corresponds to a different optimum transfer function and it is possible 
to use this transfer function as a guide in designing ttie guidance system. 
The value of this curve is that it can he used as a standard with which 
to compare the perfomnance of any system. 

Also indicated in figure 4 is the operating range of noise magnitudes 
corresponding to the class of targets and aspects previously mentioned. 

At first it might appear that each noise level would require a change in 
the filter or system so as to maintain optimum performance. However, as 
indicated hy curve B, if a system is optimized only at the mid-value of 
the noise range shown here as the design value, the performance of this 
system for other noise levels will deviate to a negligible extent from the 
optimum over the range of interest. Thus, to obtain near optimum perform- 
ance over the range of noise magnitudes likely to he encountered it is 
necessary to optimize only for the design value of noise. This is a very 
fortunate fact, and one of obvious practical importance. In a similar 
manner it can he shown that a change in the break point of the noise 
spectrum would not greatly affect the mini mum error curve shown in figiire 4 
as long as the noise spectrum is essentially flat over the handwidth of the 
optimum transfer function of figure 2. 

It is interesting to compare these results with the performance 
obtained hy disregarding noise theory in the design, that is, hy designing 
the system for the fastest possible response. Two examples are shown. 

Curve C illustrates the performance against noise which is obtained for a 
system with perfect response characteristics (unity transfer function) . 
Curve D illustrates the performance against noise for a more realistic 
system given in reference 11 (variable-incidence missile) , For this case 
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the Important satairating elements were simulated, and the system was 
designed for the fastest possible response to a step in the beam of 
100 feet. The oomparison with the optimum of the curves for both of these 
example systems shows that the perfo3rmance is significantly poorer than 
the optimimi performance indicated by theory. 


Application of Wiener Theory to the 
Beam-Rider Guidance System 


To achieve the optimum results indicated by the Wiener theory it is 
necessary to design the guidance system with frequency characteristics 
approximating the optimum specified by figure 2. The difficulties in 
accomplishing this as well as the modifications which are required are 
discussed in the following sections . 

Limitations in the application of the Wiener theory .- There are 
several restrictions on the possible forms of the transfer function which 
can be achieved in practice . The foremost of these is manifested by a 
consideration of the accelerations required of the missile. For the 
optimum system, Yq, the mean square of the required acceleration is 

expressed by 

O 00 poo 

V “ / |^w^<l>T(w)da) + / |Yo(iaj) [\4^jj(tA})du3 (9) 

-CO -00 


It is easy to see that for the form of noise spectrum given by equaticoi (l) 
and the optimum transfer function of equation (^) the spectral density of 
the acceleration required by the noise increases with frequency at the high 
frequencies, giving rise to an infinite cal led- for acceleration. Obviously 
finite values of acceleration can be obtained only if the transfer func- 
tion Yq falls off at high frequencies like l/to^ or some greater power. 

To the control designer, this restriction is perhaps more readily 
interpreted in terms of a control-deflection restriction. Assuming, for 
the sake of the argument, that the aerodynamic transfer function B/y^j 
were everywhere linear (which, of course, it is not, since practically, 
control deflections cannot increase without limit) , the analog of equa- 
tion ( 9 ) would become 


= 



+ 



= r |Yo(1w)|^ $^(w)dw + r jYo(lu3)I‘ 

'*^00 M '-Lon 




'M 


Ojj(u))dcjj (10) 
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It can te seen from equations presented in later sections that B/yw; 
approaches a constant multiple of at high, freq'uencies. Thus, 

conclusion can he drawn from equation (lO) that finite values for the 
control deflection will he obtained for any practical system which may he 
designed only if the system transfer function differs from the optimum at 
high frequencies. 

At first it would appear that any departure of the system transfer 
ftmction from the optinnim Yq 'w^uld result in increased rms error per- 
formance. However, an examination of the error spectrum indicates that, 
since the power contained in frequencies above 2 cps is small, attenviating 
the transfer function above this frequency will not affect the error 
appreciably. On the other hand, the control motion spectrum of figure 5 
shows that attenuation in the transfer ftinction above 2 cps will greatly 
reduce the required control iruDtion since the major portion of control 
motion power is concentrated in the higher frequencies. From these con- 
siderations it is obviously desirable to attenuate the transfer function 
Yq for all frequencies above the bandwidth of appreciable error spectrum. 

The attenuation might be supplied in a variety of ways. However, 
if it is assumed that an ideal filter could be Inserted in the sjrstem so 
as to supply infinite attenuation for a1 1 frequencies above 2. cps without 
affecting the response to frequencies below 2 cps, it is still found that 
for missiles similar to the one studied here deflections' required for this 
frequency band are larger than are available. Specifically, for the 
missile used herein the control deflection reqiiired is about 15 ° rms; 
since the control deflection is physically limited to about l4° (ref. 11) , 
the called-for control motions are larger than are available. Thus it 
is impossible for the actual control motion to follow the required pro- 
gram and the optimum missile motions cannot be achieved. Actually the 
value of 15 ° is quite optimistic since "the ideal filter does not exist 
and cannot even be approached very closely. Consequently, the slmul-bane- 
o\is achievement of the optimum ei*ror spectrum and called-for control 
motions which do not exceed the limited values are conflicting require- 
ments. The importance of this conflict and the solution of the problem 
posed by it are considered in the following sections. 

OptlTmnn filtering with limiting .- Because of the above restrictions 
in the application of the Wiener ■fcheoary, it is generally necessary to make 
certain modifications in the practical case. In "what follows, the effects 
of limiting and the consequent modifications required sire discussed. 

Effects of limiting: There are many limiting- type nonlinearities 

which may exist in a guidance system. Those of the most importance in 
the beam-rider systCTi are (a) control motion limiting due to mechanical 
limitations, (b) ra-fce of control motion limiting due to restrictions on 
servo capability, and (c) radar receiver voltage limi'ting due to cir- 
cuit restrictions. Typical values which are -used in these studies are 
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0.25 radian, 5 radians/sec, and 100 volts, respectively. The effects of 
these limits can be severe and need to be considered in the filtering 
problem. 

By means of Booton's recent theory 12) it is possible to evadu- 

ate the effects of nonlinearities on the over- all system performance. As 
applied to the beam-rider system, the theory shovs that it is possible to 
approximate each of the saturable elements by a simple gain which can be 
suitably chosen so that for increased rms inputs to the saturable element 
this gain is reduced. Analysis shows that in temis of the over-all sys- 
tem, characteristics the effects of this gain reduction are: (l) to reduce 

the open- loop system gain or C2) to increase the time constants in the 
open-loop transfer function. In general, both effects result in a reduc- 
tion in bandwidth of the closed- loop frequency response. In this respect, 
then, limiting effectively results in additional filtering which, of 
course, extends further into the low-frequency ranges sis the rms inputs 
increase. Thus the system response can become so slow as to result in a 
large Increase in error due to the missile’s inability to follow the t6ir- 
get motions. Furthermore, inspection of figiare 2(b) shows that the Wiener 
transfer function is conditionally stable j for large rms inputs the Wiener 
transfer function may become unstable. That this does occur in examples 
to be discussed later has been demonstrated during REAC studies. Conse- 
quently the large gain reductions and phase lags introduced by limiting 
cannot be tolerated. 

Additional filtering: In general, the energy, of the saturating quan- 

tities contained in the higher frequencies is much larger than that in the 
lower frequencies. For this reason, purposely introducing filtering into 
the system so as to attenuate the response at the higher frequencies can 
be effective in reducing the effects of limiting. This can be accomplished 
in either of two ways: by inserting additional networks in the guidance 

system, or by changing those inherent time constants which exist in any 
actual system and which are not reqiiired for the optimum filtering. 

Because of the complexity of the beam-rider guidance system there are a 
great number of places where the additional filtering can be introduced. 
These possibilities are discussed later. 

The selection of the desired frequency characteristics of the added 
filtering is more difficult than the selection of filtering location. In 
general, the choice of the desired amovint will involve a compromise between 
two extremes: (l) The first extreme is one in which the added filtering is 

too severe. As the added filtering is extended into the lower frequencies, 
limiting effects are reduced but, as previously mentioned, the results of 
limiting can only be eliminated at the expense of increased error. (2) The 
other extreme is one in which the added filtering is ins-ufficient to avoid 
serious limiting effects. This extreme is similar to the previous one 
inasmuch as the effects of limiting are also equivalent to additional fil- 
tering. For serious limiting, this equivalent filtering is too severe a nd 
reSTilts in Increased error and possibly instability. 
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It is apparent that the effects of limiting and of added filtering 
hy means of networks are closely related, in that "both are equivalent to 
the introduction of additional filtering into the system. Extremes of 
either type reduce the performance of the system siifflciently that the 
ei*ror increases due to the inability of the missile to follow the target. 


Design of the Guidance System 


In applying the modifications discussed in the preceding paragraphs 
to the design of the guidance system it is fovind convenient to consider 
two separate stages: The first stage involves the design of the system 

on a linear basis to approximate the optimum Wiener transfer function, 
while the second involves supplying the additional filtering terms to 
minimize the limiting effects. The optimum Wiener transfer function may 
be designed into the missile- control system, or into the tracking radar, 
or may be apportioned between the two, with the additional filtering in 
either place. Thus, there are a number of possibilities. Obviously, for 
linear systems a.1 1 approaches could be designed so as to produce identical 
results. In general, this does not hold for nonlinear systems , that is, 
systems in which limiting occurs, so that in the practical case it is 
necessary to consider the effect of the filtering location on performance. 
The relative merits of placing the optimum Wiener filtering first in the 
missile- control system and second in the tracking radar are discussed in 
the next two sections. 

Miss lie -control system designed for optimum Wiener filtering ,- For 
the case of missile-control-system filtering it is necessary that the 
tracking radar be designed fast enough to follow the target motion and 
the noise, since the missile-control system performs the optimum filter- 
ing operation. The first step in the design of the missile-control system 
is the synthesis of a linear system which matches the optimum Wiener 
transfer function. It has been shown that it is only necessary to approx- 
imate this transfer function over the frequency range of appreciable error. 
This can be accomplished in many ways, but to achieve this design by a 
system of conventional fonu, the design is most easily made by cut-and-try 
proced^Ires. In this proceduie it is desirable to select the aerodynamics 
as a starting point since the design of the airframe is relatively inflexi- 
ble compared to the design of the control system. The characteristics of 
a typical variable-incidence missile were chosen from reference 11. The 
transfer functions in the yaw plane for this missile are given by the fol- 
lowing equations : 


TjjS 4- 1 

T^sCTa^s^ + a^g^TaS + l) 


t 
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ZM = Tb^as + 2g-hThS + 1 

B Ts^s^CTa^s^ + 2^g^TaS + l) 

Table I summarizes the veilues of the parameters for this missile for a 
given operating condition which is held fixed in this report.- For this 
missile it is possible to synthesize a linear system to approximate 
closely the optimum transfer function. One possible system is Illustrated 
in figure 6 and the corresponding over-all open- loop transfer function is 
given by 


^ kgkaV (Tips + 1) (T4S + 1) (Tb^s^ + 2^■^^TbS + l) 

^ Td •+ kska s2(Tiis + 1) (l^s + l)(as® + bs^ + cs + l) 

The derivation of this transfer function is given in Appendix B. It 
should be noted that the transfer function of equation (ll) differs from 
that of equation (5); however, by choosing the parameters as given in 
Appendix B and as tabulated in table II, column CD , the two transfer 
functions can be closely matched over the frequency range of interest as 
illustrated by the dotted curves of figures 2 and 3 « It might be noted 
that the only restriction on the parameter T 5 in order to keep the two 
transfer functions closely matched is that it be small. Thus T 5 may 
be varied somewhat without much penalty. 

The particixlar design discussed abovei has heen chosen only for optimum 
noise performance so that other considerations important In an over-all 
design might dictate certain modifications. These possibilities are dis- 
cussed later. It is also apparent that since the design of the above 
system has been based on asstimptions of linearity the performance Indicated 
above cannot be achieved in practice due to limiting effects. As indicated 
earlier these limiting effects may result in Instahlllty. Hence modifica- 
tions are required. These modifications in the two cases of additional 
filtering, first in the tracking radar and second in the missile-control 
system, are discussed in the following paragraphs. 

Additional filtering in tracking radar: The simpler means of intro- 

ducing additional filtering into the guidance system in order to reduce 
limiting is to place It in the tracking radar since the added filtering 
and optimum Wiener filtering are achieved separately hy the tracking radar 
and missile- control system, respectively. This separation has the advan- 
tage of allowing the additional filtering to be altered without affecting 
the optimum filter design of the missile-control system. 

The additional filtering in the tracking radar may take innumerable 
forms. For this study the closed- loop transfer function of the tracking 
radar was assumed to be of the following form: 
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This transfer function -was chosen to provide "botli simplicity and a zero- 
velocity error system and it can be approximated by properly shaping the 
networks in the tracking radar. With this system, the amotint of filtering 
can conveniently be varied by a change in the time constant tx. At the 
one extreme for which tx is very small the filtering provided by the 
tracking radar is negligible, thus resulting in an unstable missile- control 
system because of limiting effects. As the time constant of the radar is 
increased, however, more fllteri33g is provided and limiting effects are 
reduced. By varying tx figure 7 was obtained; it shows the variation 
of rms error with the natural frequency, 1/2 jjtx, of the tracking radar. 

For comparison, the mini mum error obtained from Wiener theory is also 
shown. As would be expected from previous discussion, the optimum opera- 
ting point occurs at a tracking radar frequency which is within the pass- 
band of the missile-control system. At this frequency limiting effects 
still exist but are not too serious. The increase in error above the 
Wiener theory result can be attributed to the additional filtering intro- 
duced in the tracking radar and to the limiting effects in the missile- 
control system. The rms error performance of the system defined by the 
optimum operating point is given by curve E in figure 8 against a variable 
noise level. 

Additional filtering in missile-control system: An alternative place 

to Introduce additional filtering into the guidance is the missile-control 
system. This could be introduced by an additional network in the radar 
receiver. However, since it has been pointed out that the choice of time 
constant T^ was somewhat arbitrary, it is possible to increase T^ to 
provide this additional filtering. Thus limiting effects can be reduced. 

The introduction of added filtering into the missile-control system 
has the disadvantage that the added filtering also affects the system 
stability since an increase in a time lag is destabilizing. Hence, added 
stability is required. This could be provided by altering any of the 
basic parameters which are responsible for the conditional stability 
characteristic of the Wiener transfer function. To preserve the low- 
frequency characteristics of the optimum system as few changes as possible 
are desired. The parameter Txx is a convenient one with which to intro- 
duce this stability. 

The effects of the above changes were investigated by menTtH of analog 
simulation. Since in the case of missile- control- system filtering the 
tracking radar is not required to filter, it should be designed to respond 
quickly. For this purpose the radar transfer function of the form of 
equation (12) waa utilized by choosing the constants so as not appreciably 
to alter the input spectra. This was acccmiplisbed by the choice of a 
tracking radar natural frequency of 6 cps. The parameters in the missile- 
control system were altered according 'to the above discussion. Q^ical 
results obtained for the systaa optimized only at the design value of 
noise are illustrated by curve F in figure 8. The optimum Wiener perfoimi- 
ance is repeated here. 3ki general, results similar to curve F can be 
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obtained by several different combinations of parameters. One of these 
combinations requiring the fewest modifications from the Wiener system is 
given by the parameters in column @ of table II. It can be seen that 
the additional filtering is supplied by an increased radar receiver net- 
work time constant Ts while the added stability is introduced by a 
decrease in the time constant Tix. In this case the Increase in error 
above the optimum is due to the added filtering introduced by the network 
and limiting, and to the slight alteration in the Wiener system time con- 
stant Til. 

Tracking radar designed for optimum Wiener filtering .- In this section 
are considered the results ohtainahle when a.i i of the optimum filtering is 
located in the tracking radar. Here the desired optimum tracking radar 
transfer functions are given by equations (4) and ( 5 ) , and ideally the 
missile-control system should have a transfer function of unity. However, 
for many resLSons, principally those arising from, nonlinear effects, the 
latter may he expected to depart considerably from the ideal. Consequently 
it is desirable to design the missile-control system for the fastest pos- 
sible transient response within the limitations of these nonlinearities. 

Two such mis Bile -control systems were considered. The first was a con- 
ventional system chosen from reference 11 in which the response time to 
a step of 100 feet in the beam was minimized. The second was a system 
optimized for TnInimuTn response time to a small enoi:tgh step so that lin- 
earity was not exceeded. The control- system parameters for these systems 
are given in columns and ® , respectively, in table H. 

The rms error performance obtainable for tracking-radar filtering 
is summarized in figure 8 by curve G. It was found that both the 
missile-control systems gave essentially the same results. In this case 
the Increase in error above the Wiener optimum is due to slowness of the 
missile -control system and its failure to follow the beam perfectly 
because of limiting effects. 

It might appear that these limiting effects could be reduced by addi- 
tional filtering in the tracking radar. This possibility was es^lored by 
the addition of a simple first-order filter to the optimum open- loop trans- 
fer function given in equation ( 5 ) . By increasing the time constant of 
the added filter the tracking -radar performance is deteriorated frcan the 
theoretical optimum while the error of the missile-control system is 
decreased because of a reduction of limiting effects. Resialts of these 
combined effects are shown in figure 9 from which It can be seen that the 
added filtering results in progressively poorer over-all performance. Ihus 
the beneficial effects of decreased limiting are overbalanced by the detri- 
mental effects of altering the tracking -radar transfer function fran the 
optimum. 



NACA RM A55E11 


IT 


CcanpariBons and Other Considerations 


The results of various filtering arrangements are ccmpared in fig- 
ure 8 on the basis of the ims error. If it were not for the nonlinearities 
the performance for both missile -control system filtering and tracking- 
radar filtering coiild be made identical with the performance of the optimum 
Wiener system. However, it is interesting to find that even in the non- 
linear case ai 1 arrangements can resxilt in ccai 5 »ar*able error performance. 
Some advantage in tracking- radar filtering is apparent from the figure. 

For comparison, again, the performance of a typical missile- control system 
which was optimized for fast response in the absence of noise is repeated. 

Up to this point only the error performance has been considered. 

Other factors which are of importance in the over-all evaluation of a guid- 
ance system will now be discussed. One such faetor is the servo energy 
required to achieve a given error performance. For a servo system al r eady 
designed to meet the maximum expected hinge moment, the servo power is 
proportional to the time average of the sum of the absolute displacements 
of 5 between values at which 6 changes sign. Thus the average servo 
power for a time t is 


Z ££> 

average servo power — 

Evaluation of the servo power for the systems discussed above has shown 
that both methods of missile- control-system filtering require about 
22 percent more power than does tracking-radar filtering. Thus in an 
over- all evaluation based both on rms error performance and servo energy 
requirements it is apparent that tracking radar filtering is slightly 
superior . 

There may be still other factors of importance in guidance system 
design even within the framework of the assumptions discussed earlier. 
Usually these requirements are related, to the specific design objectives 
of the system and may dictate certain modifications such as the choice of' 
filtering location or ailteratlons in certain individual transfer functions 
of the system. For example, for short-range missiles in which launching 
eirors are not prevented frcan building up, capture of the beam in Tn-f-n-TTmim 
time may be important enough that tracking -radar filtering would be pref- 
erable to missile-control-system filtering because of the fast response 
which can be designed into the missile -control system. In other cases 
it may be desirable to alter the design of the system somewhat for any of 
a number of reasons. Requirements of simplicity on certain parts of the 
system or the necessity of using certain fixed and unalterable elements in 
the design are examples. Another possibility is that for flight conditions 
inT7hich serious atmospheric turbulence exists it might be necessary to 
make certain alterations to minimize the response to gusts.. An investiga- 
tion of such factors is beyond the scope of this report. 
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CONCLUDING REMARKS 

This study has considered the application of Wiener filter theory to 
the optimization of a beam-rider system operating in the presence of noise. 
The Wiener theory hae been used to establish both the theoretical lower 
limit of error and the desired transfer-function characteristics. Although 
these transfer-function characteristics might be ejqpected to vary with 
noise magnitiide, it has been found that when the transfer function is cho- 
sen to be optimum for only a Md-value of noise the performance for other 
noise levels deviates to a negligible extent frcan. the optimum. 

In general, linear beam- rider systems can be synthesized to produce 
theoretical optimum performance, but the practical achievement of these 
results is restricted by limiting of control deflections, rate of control 
deflections, and radar receiver voltages. With suitable modifications, 
however, the theory can be useful as a guide in system design. Optimiza- 
tion in the actual nonlinear case is shown to involve two considerations : 
limiting and additional filtering. It is shown that the important types 
of limiting tend to result in system Instability. However by appropriate 
placement of additional network filtering it is possible to minimize these 
limiting effects without serious deterioration of the error performance. 

The design of a guidance system is most conveniently accomplished in 
two stages: The first consists of designing the system on a linear basis 

to approximate the Wiener transfer function, and the second of supplying 
the additional filtering terms to minimize the limiting effects. The 
application of the Wiener theory and the modifications required to arrive 
at an optimum system design have been illustrated in this report by con- 
sidering systems in which the optimum Wiener filtering is designed into 
either the miss lie -control system or the tracking radar, and additional 
network filtering supplied in either place. It was found that comparable 
error performance can be achieved by any of the methods. Consideration 
of both rms error and servo energy requirements for the cases studied 
indicates that optimum filtering in the tracking radar is slightly superior 
to that in the missile -control system. 

In modifying the results of the Wiener theory in order to minimize 
the effects of limiting, the best results were, in general, obtained when 
the additional filtering was added in such a way as to tend to keep the 
operation of the system in the linear range. This has suggested the pos- 
sibility of seeking an optimum solution based on the ' stipulation that the 
filtering should restrict the operation of the system to within its linear 
range. In reference 13 the results of such an analysis are presented and 
compared to those of the present report. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 11, I 955 
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APFEHIHX A 

DETEEMUtATION QP THE OPTIMUM WIENER EILIER 


Wiener’s solution for the optimum transfer fimction Yq can he 
e3^)ressed hy the following equation: 


Yo(ico) 


. J-... r"e-io.t r 

2n$‘^(co)Jo '^oa ® (c.) 


do, dt 


(Al) 


where represents the target motion spectral density, and the quanti- 

ties and are defined as the factors of a certain function $ with 
poles and zeroes located in the upper and lower half -plane a, respectively; 
the factor ® is defined as equal to <5^ + so that 

$ = + •Bjj' — (A2) 

For the case to he considered here, the target motion and noise are 
defined hy 


itco^(io^+ k^) 


®N = N 


(A3) 

(Alt) 


It will he noted that the noise spectrum has heen approximated here hy a 
constant in order to reduce the complexity of the calculations. This 
approximation is valid because the more exact form of the noise spectrum 
(eq. (l) ) is essentially flat over the handwidth of the optimum transfer 
function. Use of the more exact form of the noise spectrum affects only 
the response at the hi^ frequencies which are heyond the range of inter- 
est. It should also he pointed out that the use of equation (A3), as 
such, leads to certain mathematical difficulties in evaluating the right- 
hand side of equation (Al) , hecatise the theory requires that the poles 
of ®ip not he located on the real axis. Rigorously, to avoid these 
difficulties, it is necessary to modify the target motion spectral density 
to the follo^ring: 


ffirn = 

^ jc(ai2 + Tll^) (*1^ + Tl2^)(^ + 

where T^x snd r\ 2 . any small real numhers. Thus the solution of equa- 
tion (Al) is a f\mction of ioj, Tix/ ami t\z- desired sinswer is then 

obtained as the limit of Yo(lcd,'nx^‘na) as t\x snd Tja approach zero. 
However, it can he shown that the same answer can he obtained more sin^jly 
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by taking tii and Tja equal to zero as vould Tae obtained frcm the rigor- 
ous process described above. It will be convenient to retain the ti’s 
for a few lines; subsequently they will be dropped. According to equa- 
tion (A2) then, 

4. = H g° + + (ka==M) 

(a + iTii)(a + iTja) (a - ii^i) (a - ±r[z)(a.^ + k^) 

Ca6) 

(g - ai) (g - Og) (a - 03 ) (<x + ai) (a, + oog) (g + oq) / ._v 

(a + iTji) (g + 1 t] 2 ) (g - lTii)(g - iTi 2 )(g + ik) (g - ik) 


where gx^ cx^f an d oq represent roots in the upper half -plane. It should 
be noted that none of these roots are real. From equation (A2) also, it 
follows that 

«+(g) = N “ ^x)(a‘ - a2)(g - 03 ) ^ ^ p(a.) 

(g - iTii)(g - 1^2) (a - ik) (a - iTii)(g - i'q2) ■" 

(A8) 


«Ji(a) = 


®(g) = 


$"(g) 


_ (g + gx) (g ocffi ) (g + g<3 ) _ 


“P(-g) 


(g + iTix) (ct + 11)2) (g + ik) ~ (g + Iqx) (00 + 1112) (a + Ik) 


(A9) 


According to the definition given above 

p(g) = (g - ax)(g - 02) (a - 03) 
s g® + b2U^ + bxg + bo 

and 


(AlO) 


-p(-a) = g® - b 2 cc 2 + bxtt - bo 

where 

b2 = “(ax + Os + o£e) 
hi = aict2 + ai 03 + 0203 
bo = -aio;2°<3 


(All) 
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Tile Inner integral of equation (Al) is 




eig.t 

- ik) (a, + ai) (a + 02 )Ca + oq) 


dot 



(A13) 


In the above expression it should be noted that it has been assumed tii 
and are zero in order to simplify the following expressions . Now if 
a is considered to be a complex variable, the above integration is 
equivalent to integration over the contour shown below. 



Here the only two poles involved are a second- order pole at the origin 
(actually at iTii and ir^a in the rigorous case) and a first-order pole 
at a = ik since, as indicated above, -ai, -oaj and -oq lie in the 
lower half -plane. Hie pertinent residues can be found in the usual man- 
ner as shown below: 


Res(o) 


Htti 
a o 




t _ ikbi + bp 

Icbo k^o^ 


(Al4) 


Res(ik) = T -Ttti [ (o, - ik)f (a) ] 
a ik 

e~kt 

-k^(ik + ai) (ik + Oa) (ii*^ + °e) -lJ^[-p(-ik) ] 


(A15). 
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Hence the integral in equation (AI 3 ) hecomes 


Ii = 2ka^i 
= 2k&^l 


{- 


ikbi + ho 
k^bo^ 


71 + Sfe + 


kbo -k2[ 
1 


e-kt ^ 
-p(-ik)i; 


(A 16 ) 


■where the definitions of y ^ and 7 ^ are appaxent. 

The second Integration in equation (Al) is merely a Foxirier transform 
and is found as follo-ws: 


/ OO 

e-i“'blidt 


= 2ka^i 


r 


7,e“i‘^'f^ + rS- + 7^e"(^ + 

^ J^Dn ^ 


dt (AI 7 ) 


Thus, 


I 2 = 


2ka2i‘ 


LitJ kbo(io))2 k 


+ icjJ 


= 2ka^i 


(7i + 7s) (±^)^ + (k7i + 

(i(jk>)^(k + ico) 


(AI 8 ) 


It should be no-ted that because and t\z 'were assumed to be zero for 
simplicity, questions concenaing the existence of the above in’tegral 
arise. Ho'wever if the analysis is made 'without this assun^itlon, it can 
be shown that the integral in equation (AI 7 ) does exist, and that the 
limit of this integral as t)i and go to zero becomes precisely equa- 
tion (A18) . The coefficients in this equation can be simplified to a 
more useful form by the following development. From the definitions 
given earlier in equations (A6) through (A9) , 

a® + k^a'*' + = [p(a)][-p(-a)] 


which at a = Ik becomes 

^ = -bo^ = [p(ik)][-p(-ik)] 
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When this relation is substituted in the definition of 7 ^ given in 
equation (AI 6 ) , 

1 _ 

~ -l?[-p(-lt)] (-k=) (-bo^) 


-ll£? - bak^ + ibik + bo 
" ^b^ 

Thus the coefficients in equation (AI 8 ) simplify to 

ik + ba 


7i + 7-5 = “ 


130-= 


The transfer function Yq cen now be found from the preceding 
equations as 


Yo(ico) = 


. ^ 1 


2jt*+ 


(- * + ^X' ^ OG ^ 0 


(A 19 ) 


vhich, in terms of the conventional complex frequency s = iu, reduces to 
the following alternative forms: 


Yo(s) = 


ik + b 


"bo 


■> • (- =&> 




+ 1 


Yo(s) 


Yo(s) = 


(» • i)(# • iX» • "■) 


+ 1 




■5^ 


s - 




(A20) 


o ^o “o 


+ 2ggTt;^a + 1 

(TpS + 1) (Ty^B^ + ^iyTyS + 1) 
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Qh 


where 




Ih + ha 


ho 




a 


bo 


qiZ _ _ — 1 

•*■7 oeciQ 






tti 


In terms of the equivalent open-loop transfer function^ \^q, of a unity 
feedback system, 

^ 1 -4o(sy 

Ks ^ 0 


Ik 


_ Ig"°" - ggg'JaB 1 , 

■ ^ s2(t^8 + 1) 


where 


^ ” tS 

The coefficients of the transfer function of the optimum system have 
been evaluated for the following values of target motion and noise: 

k = 0.4 

T = 5 sec 

a = 1 g 

N = 15 ft^/radlan/sec 

Evaluation of the numerator of equation (a 6) gives 


+ k^a** + 


a® + 0.l6 a'*' + 8.8 

(a - o,i)(a - 02 )(a - oq) (a, + ax) (a + Og) (a + oa) 
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where 


ai = i 1.456 

Qte = 1.427 exp(i O. 543 ) 

as = *1.427 exp(-i O. 543 ) 

■which. Eire e-l ~l loca'ted in "bhe upper half-plane. These roots are then 
sufficient to djetennine the optimum transfer function. The constants in 
equation (A12) give 


ha - -± 2.94 
hi = -4.201 
ho = i 2.97 

Then evaluating the parame’ters helow equation (A20) and in equation (A2l) 
gives the following values: 


Parame'ter 

z 


Value 

0.925 



.700 

.519 

2.50 

7.42 


Frcm these parame-bers the optimum transfer function of equation (A20) 
heccm.es 


Y ( 0.855 s2 + 1.4l s + 1) 

( 0.687 B + 1 ) ( 0.490 S^ + 0.727 s + 1) 


(A22) 


and from equation (A2l) the equivalent open- loop transfer function hecomes 

,,(s) . 7.42 ° 

s^( 2.50 s + 1) 


(A23) 
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APFEHDIX B 

MISSIIE-COMTROL-SYSTEM APPROXIMA!nO!f TO THE OPTIMUM 
TRANSFER FUNCTION 


An optimum linear design Is illustrated in figure 6 , and in the fol' 
lowing the system equations are derived. Frcan figure 6 , 



ksCT^s+l) (T 3 S+I) 

Tds(TxS+l) (Ta2s2+2^a.TaS+l) (TaS+D+ksl^sS (l^js+l) (T 2 S+I) 


(Bl) 


To simplify equation (Bl) it is convenient to make T 3 = Tj^. This choice 
is not essential hut its use leads to more easily handled equations. In 
certain cases where gust disturhances are serious it may be more desirable 
to choose T 3 small. With the former assunq)tion, however. 



ks(TMS+l) 

s [ (TdTa^Ti) s3+(TdTa2+Td2^aTaTi) s2+(TdTi+Td2^aTa+ksl£3T2 ) s+(Td+ksh3 ) ] 


kg Tm»+1 

Td+kgka s (as^+bs^+cs+l) 


(B2) 


where 


•'1 

Td + kska 

Td(y + 

Td '+ kgks I 


25g^TaTd + TiTd + kgksTs 
° ^ Td + kgka 

Thus the entire open-loop transfer fiinction can he written as 


(B3) 


_ ^ = kgkgV (T 4 S + 1)(Tiqs 4 - l)(Tb^s^ + 2g-(3Tb8 + l) 

^ I Td + kgka s^(TgS + l) (TiiS + l) (as® + bs^ + cs + l) 


It is desirable to choose the parameters in this equation so as to match 
the optimum transfer function given in equation (A21) . Since the system 
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equation (Bh) is of much higher order than the optimum, it is apparent 
that the matching cannot he acccanplished perfectly. However, as shown 
in the text it is only necessary to app3X)ximate the optimum transfer func- 
tion over a limited frequency hand, Ihis can he accomplished hy choosing 
certain terms in (b 4) to correspond to teims in (A21) . The remaining 
terms should then he chosen to have negligible effect on the optimum error 
spectrum. 

There are many ways the design can he accomplished. No attempt will 
he made to investigate all the possihillties. One design, however, is 
based on the following correspondences: 

T4 <—> Tb, 

^10 ^ ^ 

Til Tx 


The remaining terms can he chosen in any manner as long as their effect 
is small over the optimum error spectrum. One of the many possible 
choices is to factor the cubic to approximate the aerodynamic term in the 
numerator as follows: 


as® + hs 2 + cs + 1 = (Tb®s® + 2 ^i 3 Tbs + l) (T-^s + l) (B 5 ) 


where Tj^ is arhitraiily chosen to he small. This factoring is not 
essential to the design hut it leads to simple equations for the control- 
system parameters. For example, for given aerodynamics and, choice of 
the coefficients of the left side of (B 5 ) are determined hy 


a = T^Tb^ 

h = TipCi^Tb + Tb® ► 
c = Tl + a^T^Tb 


(b6) 


Solving equation (B3) tt^n for the control- system parameters gives 

T^Tb^Ta ^ ^ 

" 2CbTbTiTa + Tb^(Ta - ~ ^ 








2Tb" 


^^2 _ ^2 f^h^ " Sa^a) + Ti + 2 ^ijjTb 


(B7) 
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Then the value of kgks is determined froni the desired gain; 


kaks = k^ 


It will be noted that since only kaks and keka are specified, one of 
the three may be chosen arbitrarily* The numerical values of the above 
parameters have been calculated and are listed in column @ , table II. 
It should be pointed out that this design is only one of many possible 
designs. In any particular case cearbain modifications may be necessary 
as discussed in the section on ccmparisons and other considerations . 
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TABIE I.- SUMMART OF AERODYNAMIC PARAMETERS FOR EXAMPLE MISSILE 


U) 

o 



TABLE II.- SUMMARY OF CONTROL-SYSTEM PARAMETERS 


Parameter 




Linear 

mis sile- control 
system optimum 
with noise 


0.25 

.912 

.01593 

.912 

2.5 

.0it62 

.01603 

.0122 

.81^6 


Nonlinear 
ml s s lie - c ontrol 
system optimum 
with noise 


0.25 

.912 

.21^9 

.912 

.80 

.01603 

41^.2 

.0122 

.846 


© 


Nonlinear 
missile -control 
system optimum 
for 100-foot step y 3 


1.0 

.32 

.0559 


.0358 

.025 

44.2 

.0721 

.846 


Linear 

mi s sile - c ontrol 
system optimum 
for step ■ yg 


.0358 

.025 
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Frequency, cps 

(a) Closed -loop transfer function, Y*. 

Figure 2.- Optimum transfer functions. 


Phase angle, degrees 
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Frequency, cps 

(b) Open-loop transfer function, /io* 


Figure 2.- Concluded. 


Phase angle, degrees 


34 


NACA RM A55E11 


(A 

=] 

o 


U> 

+“ 

c 

0) 

c 

o 

Q. 

E 

o 

o 


o 

L. 

0) 



OT 

c 

a> 

•o 


o 


o 

<u 

Q. 

0> 


O 


o> 

Q 



Frequency, cps 


Figure 3,-. Spectral density of error components. 




RMS error, JW ^ feet 


NACA RM A55E11 


3 : 






T.V.V.'.WV.V.V 
!¥•••■• w^y* • • • • * 

>**mA.O****>» 


a • • a • • ■«•>•••• ca>Aa a ••••••««« 




« |-V*V.V.JV,V»JiV.W 

• •«•■»•••••■•■ •••■*• •%%*•*•*•*•* »*aVa%%*^%VjiVa 

aaaaa»aaaaaaaaaaaaaaaaaaaaaaaaaaa»aaaaaaaa 
•aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaai 


Optimized for no noise 

System with limiters 

System with unity 
transfer function 






'iaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

laaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaiaaaaai 

ia»aaasaaaaasaaaaaaaaa»aaaaaaaaa»aaaaaaaaaaaa»a 

aaa»aaaaaaa»aaaaaaaaaaa»aaaaaaaaaaaaaaaaaaa»aai 


'aaaaaaaa»aaaiaa*aaaaaaaa»aaaaaaaaaaaaaaaaaaaaa 
aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaai 
a\Va\*aVaVaVa\'a^a*'a*a*«*aV\* *^***\**J* ***************** 
'aaaaaaaaaaaaaaavaaaaalaaaaaaaaaaaaaaaaaaaaaava 

la aaa aa aa*a aaa at a aaaaa a aaa aa»aaaaaaaaa*a*aaa»aa 


■■aataaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 
■ aaa aaaaaBaaaaaaaaaaaaaaaaaaaaaaaaa*aaaaaaaaaaa 

/Ivivivlw 




^Operating range 


• aaaaaaaaaa a a a a a a a a a a a a a a a a a a a a a a a a a a a < 


aaaaaatBaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa aajTa a a aa 
• aaaaBaaaaaaaaaaataaaaaaaaaaaaaaaaaa a^Taa a a a a a a 



Optimized for design 
value, Wiener theory 


Minimum error, 
Wiener theory 


Sp^^^Design value 




Noise spectral density magnitude, N, ft^/nadion/sec 


Figijre U.- Optimvnn performance 137 Wiener tiieory. 


spectral density of control-motion components, lorus 
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Frequency, cps 


Figure Spectrai density of control motion components 
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Figure 



7,~ Effect of adjcLitlonal filtering in tracking radar; missile- 
control system designed for optimum Wiener filtering. 
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E. Missile-control system designed for optimum Wiener 
filtering with additional filtering in tracking radar 

F. Missile-control system designed for both optimum 
Wiener filtering and additional filtering 


G. Tracking radar designed for optimum Wiener filtering 



Figure 8.- Effect of practical limitB on minimum ein*or; 

design valiie. 


optimum for 
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ho 



Time constant of added filter, seconds 


Figure 9 .- Effect of additional flltOi*ing in tracking radar designed for 

optijarum Wiener filtering. 


NACA - tongley Field, V». 





